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ABSTRACT 

Lithium abundance derived in metal-poor main sequence stars is about three times 
lower than the value of primordial Li predicted by the standard Big Bang nucleosynthe¬ 
sis when the baryon density is taken from the CMB or the deuterium measurements. 
This disagreement is generally referred as the lithium problem. We here reconsider the 
stellar Li evolution from the pre-main sequence to the end of the main sequence phase 
by introducing the effects of convective overshooting and residual mass accretion. We 
show that ^Li could be significantly depleted by convective overshooting in the pre- 
main sequence phase and then partially restored in the stellar atmosphere by a tail of 
matter accretion which follows the Li depletion phase and that could be regulated by 
EUV photo-evaporation. By considering the conventional nuclear burning and micro¬ 
scopic diffusion along the main sequence we can reproduce the Spite plateau for stars 
with initial mass mo = 0.62 — 0.80 Mq, and the Li declining branch for lower mass 
dwarfs, e.g, mo = 0.57 — 0.60 Mq, for a wide range of metallicities (Z=0.00001 to 
Z=0.0005), starting from an initial Li abundance ^(Li) = 2.72. This environmental Li 
evolution model also offers the possibility to interpret the decrease of Li abundance in 
extremely metal-poor stars, the Li disparities in spectroscopic binaries and the low Li 
abundance in planet hosting stars. 

Key words: star: pre-main sequence - abundance: lithium. 


1 INTRODUCTION 


^Li is one of the four isotopes synthesized in the primordial 
nucleosynthesis. Its primordial abundance depends mainly 
on the baryon-to-photon ratio with only a minor sensitivity 
to the universal speed- up expansion rate, namely the num¬ 
ber of neutrino families (lFieldsll201ll : lFields et al.l2014l ). The 
universal baryon density can be obtained either from the 
acoustic oscillations of the cosmic microwave background 
(CMB) observation or independently, from the primordial 
deuterium abundance measured in un-evolved clouds of 
distant quasar spectra (lAdamI [l976l ). Observations based 
on the Wilkinson M icrowave Anisotropy Probe (WMAP; 
iKomatsu et 1201 ll) predict a primordial ^Li abundance 
A(Li) = 2.743 dCoc et al]l2012f). From the baryon density 
measured by the Planck mission dPlanck Collaboration Xv3 
l201,'j l ICoc et ^ l|2014l ) calculated the primordial value of 
■^Li/H to be 4.56 - 5.34 x 10"^° (A(Li)« 2.66 - 2.73). 


* E-mail: xtfu@sissa.it 

^ A(Li) = 12-1- log[n(Li)/n(H)] where n is number density of 
atoms and 12 is the solar hydrogen abundance. 


POP II m ain sequence (MS) stars show a constant 
^Li abundance dSpite &: Spitd[l982l ). which was interpreted 
as an evidence that the these stars carry the primordial 
^Li abundance because of their low metallicity. In the 
past three decades, observations of metal- poor main se¬ 
quen c e stars both in the M ilky Way halo l|Spite &: Spitd 
1 19821: ISb ordone et all l201f)l ). and in the globular clusters 
llLind et al.ll20od) conf i rm th at the ^Li abundance remains 
A(Li) « 2.26 ( Molard I2n08h . This abundance, which de¬ 
fines the so-called Spite plateau, is three times lower than 
the predicted primordial value. The discrepancy is the long¬ 
standing “lithium problem”. 


There are several lines of study which have been pur¬ 
sued to provide possible solutions to the problem: i) nuclear 
physics solutions which alter the reaction flow into and out of 
mass-7 dCoc et al ][2 m 3); ii) new particle p hysics where mas¬ 
sive decaying parti cles could destroy ^Li (lOlive et al.ll20l3 : 
iKaiino et al.l 1^121 ): iii) Chemical separation by magnetic 
field in the early stru cture formation that reduces t he abun¬ 
dance ratio of Li/H (iKusakabe &: Kawasakil [20151) : iv) ^Li 
depletion during main sequence evolution. It has been ar¬ 
gued that certain physical processes, which may occur as the 


© 0000 RAS 





































2 X. Fu, et al 


stars evolve on the main sequence, could cause the observed 
lithium depletion. A mong these proc e sses w e recall gravita¬ 
tional settli ng fe.s. ISa laris fc Weisd ll200lll : iRichard et al.l 
1I2OO5II and iKorn et ah ( 2006l ll or possible coupling be- 
tween internal gravity wave s and rotation-induced mixing 
llCharbonnel fc Prima.j|2005l l. 

In this paper we address the lithium problem on 
the ground of the Pre-Main Sequence (PMS) stellar 
physics, under the assumption that the standard BBN 
is correct. Hereafter when not otherwise specified for 
Li we refer to the isotope 7. The surface Li evolution 
during the PMS has been analyzed in stellar evolu- 


dD’Antona & Mazzitellil 1984: 

Soderblom & Fedele et al. 

I993I: iD’Antona & Mazzitellil 

Il994: 1 Swenson et al. 

I994I; IVentura et al.l 19981: Piau & Turck-Chiezd 20021: 


Tognelli et al.l 2012h . with the general conclusion that 


mixing efficiency plays a key role on Li depletion. Recent 
extensive studies in metal-rich young open clusters do show 
that the observed lithium abundance is modifie d during 
the PMS phase (ISomers fc Pinsonneaulti |2014|L In this 
respect, two points are worthy of consideration. On one 
side, it has been shown that overshoot at the base of the 
solar con vective envelope is more efficient than hitherto 
believed JChristensen-Dalsgaard et al.ll201ll L and could be 
even more in the envelope of a PMS star. This would favour 
a more efhcient Li depletion. On the other side, there is 
evidence that a residual mass accretion (also called late 
accretion) persists for several 10^ y ears and is observed u p 
to the early main sequence phase (iDe Marchi et ah llml). 
This may act to partiall y restore Li in metal-poor PMS 
stars. iMolaro et 'aU ll2012tl first proposed that Li abundance 
both in the Population I (POP I) and Population II (POP 
II) stars could be modified due to the combined effects 
of efficient overshoot and late mass accretion during the 
PMS evolution. In this study we use PARSEC (PAdova and 
TRieste Stellar Evolution Code; iBressan et al. 2012ll to 
exami ne the working scenario proposed by Molaro et al.l 
ll2012ll in a more systematic way. 

The structure of the paper is as follows. Section 2 de¬ 
scribes the theoretical pre-main sequence models, which in¬ 
clude envelope overshooting, residual accretion and EUV 
photo-evaporation. Section 3 describes the Li evolution in 
main sequence. Section 4 presents the results of the models 
and the comparison with the observations. An ample discus¬ 
sion and the main conclusions are drawn in Section 5. 


2 PRE-MAIN SEQUENCE LI EVOLUTION 

Pre-main sequence is the direct continuation of the proto- 
stellar phase. Initially, as the young stellar object (YSO) 
evolves along its stellar birthline, its luminosity is mainly 
supported by an accretion pr ocess strong e nough to main¬ 
tain active deuterium fusion dStahlerj ll983l L Once the ac¬ 
cretion ceases, the proto-star, surrounded only by a residual 
disk, descends along Its Hayashi line almost vertically In the 
Hertzsprung-Russell diagram (H-R diagram), undergoing a 
rapid gravitational contraction. 

In this phase evolution of PMS stars critically depends 
on their mass, and lithium can be burned at different stages 
through the reaction ^Li(p,a)'^He. Lithium is very fragile as 


the nuclear reaction rate Rnuciear of ^Li(p,Q:)'*He becomes 
efficient already at temperatures of a few million Kelvin. 
The effective Li burning temperature for PMS stars, i.e. 
that needed to consume Li in a timescale of ~ 10^ yr, is 
~ 4 X 10® K. The Rnuciear is adopted from JINA REACLIB 
database dCvburt et ahll^lOh . 

Very low mass stars with initial mass mo < O.O6M0 
never reach this temperature. More massive PMS stars ex¬ 
perience Li-burning, which initially affects the entire stellar 
structure as long as it is fully convective. Later, at the for¬ 
mation of the radiative core, the extent of Li burning can 
vary, depending on the mass of the star. Thereafter, the Li 
evolution is critically affected the temperature at the base 
of the convective envelope, and the efficiency of overshoot 
at the base of the convective envelope begins to play a sig¬ 
nificant role. 

We use the stellar evolution code PARSEC to calculate 
the PMS evolution of low mass stars with initial mass from 
0.50 Mq to 0.85 Mq, with metallicity Z = 0.0001 (~ 
\M/H] = —2.2 dex), typical of POP II stars. The adopted 
helium abund ance is Yj, = 0.249^ based on the helium en¬ 
richment law dBressan et al.l[2012l') : 

AT 

y = yj,-P-—z = 0.2485 -fl.78Z (1) 

LXZj 

where Y), = 0.2485 is the adopted primordial value and 
/ l\Z is the helium-to-metals enrichment ratio. Following 
the SBBN prediction, we set the initial lithium abundance 
A(Li)=2.72. A solar-model calibrated mixing length param¬ 
eter umlt = 1.74 is adopted. 


2.1 Convective overshooting 


Overshooting (OV) is the signature of the non-local convec¬ 
tive mixing in the star that may occur at the borde rs of 
any convectively unstable region dBressan et al.l[2014l L As 
already anticipated, efficient overshooting at the base of the 
convective envelope may significantly affect PMS surface Li 
depletion, because overshooting extends the mixed region 

into the hotte r stellar interior. _ _ 

Recently Christensen-Dalsgaard et al.l (1201 il l suggest 
that an overshooting region of size Ae « 0.3 ~ 0.5 Hp 

(where Hp is the pressure scale height), at the base of the 
convective envelope of the Sun, provides a better agreement 
with the helioseismology data. These values are higher tha n 
those of earlier studies (e.g., iMiglio fc Montalban.1 I^OSll . 
and indicate that the overshoot mixing could be much more 
efficient than previously assumed. It should be recalled that 
by applying an efficient envelope overshoot, many other ob¬ 
servational data could be better reprodu ced, such as the blu e 
loop of intermediate and mass ive stars llA longi et al.l ll99lll . 
the low-mass metal-poor stars dBressan et al. 2012ll . the car¬ 
bon stars luminosity fnnctions in the Magella nic Clouds via 
a more efhcient third dr edge-up in AGB stars (iHerwi ^ l2000l : 
iMarigo fc Girardil |2(I07 ), and the location of the RGB bump 
in old star clusters dAlongi et al.lll99ll L 

Soon after the formation of the radiative core we as¬ 
sume that the overshooting at the base of the deep convec¬ 
tive envelope is quite efficient, with Ae ~ 1.5 Hp. This 

value is intermediate between the one suggested for the 
Sun dChristensen-Dalsgaard et al.ll20U) and the larger val¬ 
ues (Ae = 2 — 4 Hp) suggested by Thng et ahl d2014h for 
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Figure 1. Kippenhahn diagrams for stars with initial metal- 
licity Z=0.0001 and different initial masses. The hatched area 
corresponds the mass size of the convective zone, within it the 
area filled with vertical lines represents the contribution of the 
envelope overshooting. On the y-axis the mass coordinate of the 
convective zone starts from 0 (at the center of the star) which 
means the star is fully convective, while 1 means the base of the 
convective zone is at the surface of the star and there is no convec¬ 
tion. As the star evolves, the base of the convective zone retreats 
towards the surface of the star, and the overshoot layer become 
shallower and shallower. 


intermediate and massive stars. We will show in Sect. [5] 
that this value is not critical, since even assuming a max¬ 
imum overshoot Ae ~ 0.7 Hp, the model can well re¬ 

produce the data. The overshooting distance is computed 
using the natural logarithm of the pressure, downward 
from the Schwarzschild border (PscJi) at the bottom of the 
convective envelope, i.e. an underlying stable layer (Pjay) 
is mixed with the envelope if In(Piay) - In(Psc^) Si Ae. 
Then, as the star moves toward the zero age main se¬ 
quence (ZAMS), we vary the overshooting efficiency pro¬ 
portionally to the mass size of the outer unstable convec¬ 
tive region (fez), un til the typical value estimated for the 
Sun Ae = 0.3 Hp lIChristensen-Dalsgaard et al.l l 2 nilh is 
reached: 


Ae = 0.3 + (1.5 - 0.3) > 1 = /e; 


( 2 ) 


This assumption produces an efficient photospheric Li de¬ 
pletion during the PMS phase, while preserving the solar 
constraints for the main sequence stars. Fig. [T] illustrates 
the temporal evolution of the size of the convective and 
overshoot regions for four values of the initial stellar mass, 
mo = 0.5, 0.6, 0.7, 0.8 Mq at metallicity Z=0.0001. The mass 
coordinate, on the Y axis, goes from 0, at the center of the 
star, to 1 at the surface. The stars are initially fully con¬ 
vective and the temperature at the base of the convective 
zone {Tbcz) is high enough to burn Li so efficiently to be¬ 
gin surface Li depletion. As the central region heats up, the 
base of the convective envelope moves outwards to the stel¬ 
lar surface, while Ti^cz decreases. When Tb^z drops below 
the lithium effective burning threshold, the photospheric Li- 
depletion ceases, leading to an almost total or only partial 
Li depletion, depending on the stellar mass and on the over¬ 
shooting efficiency. 

Figure [2] illustrates the effect of efficient overshoot- 
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Figure 2. PMS Li evolution in the standard model (left panel) 
and model with envelope overshoot (right panel) for three stel¬ 
lar models under consideration. The two horizontal lines indicate 
primordial A(Li) (dark grey dot dashed line) and Spite plateau 
A(Li) (light grey solid line), respectively. Overshooting is the solo 
parameter tested, no other mechanism is applied. 


ing as the solo mechanism applied to the PMS phase for 
mo = 0.60, 0.70, 0.80 Mq stars. In standard model without 
overshoot (left panel), A(Li) almost remains the same as the 
primordial value, unless the convective zone in the star itself 
goes deep enough to burn Li, as for mo = 0.60 Mq star in 
the left panel. Conversely, assuming efficient envelope over¬ 
shooting, Li abundance is found to be fully depleted over a 
few Myr in the same stellar models (right panel). It follows 
that if overshoot were the only process at work, a very low 
Li abundance should be measured at the surface of those 
stars. 


2.2 Late mass accretion 

Recent observations of late PMS stars - those that have 
already abandoned the stellar birth line and are joining 
the ZAMS - reveal that most of them show Ha-line emis- 
sion, which likely originates in a residual a ccretion process 
llDe Marchi et ail l2010l : ISnezzi et al.l |2012| ) . This suggests 
that the disk accretion may last much longer than previ¬ 
ously believed: at lea st tens of Myr and even up to the early 
main sequence stars dPe Marchi et aT]l201ll ). 

This late accretion (also known as residual accretion) is 
different from the one that maintains the star on the stellar 
birth line during the proto-stellar phase. The observed val¬ 
ues indicate a median accretion rate of (~ 1 O“®M 0 /yr) for 
YSO T Tauri stars (TTS) w ith a steady decline as time pro- 
ce eds dE spaillat et al.l 2013 ). Most residual d i sks p resented 
in lEspaillat et al.l (l2014h and ICallardo et al.l l| 2 nid ) are as¬ 


sociated with accretion rates at the level of ~ 10 '^Mq/ji 
- lO"®M0/yr. 

Since the accreting material keeps the initial Li abun¬ 
dance, its contribution cannot be neglected, especially if the 
photospheric value has been already depleted by another 
process. In this framework we modify our stellar evolution 
code PARSEC to account for the effect of such a residual ac¬ 
cretion during the PMS phase. 

We assume that, after the main accretion phase when 
the star leaves the stellar birth-line by consuming its inter¬ 
nal deuterium, a residual accretion keeps going on. Since our 
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Figure 3. An example of the PMS lithium evolution for a star 
with mo = O.75Af0, |M/H|=-2.2. Left panel: Evolutionary track 
at constant mass in H-R diagram. The solid line starts at the 
end of the dotted-line-indicated deuterium burning phase. At this 
stage the accretion is reduced to 2 x 10“® Mq/jv. The asterisk 
marks the end of the residual accretion. Right panel: Li evolution 
starting from an initial abundance A(Li)=2.72 dex (horizontal 
black dot dashed line). The temperatures at the center (Teen, 
dark red dot dot dashed line) and at the base of the convection 
zone (T^cz, dark red dashed line) are also shown. The accretion 
(with rate drawn by the dark blue dotted line) is terminated by 
the EUV photo-evaporation (with rate drawn by the dark green 
dot dashed line) when the latter reaches the same value as the 
former. 


PMS models are initially evolved at constant mass from a 
contracting configuration without nuclear burning, we as¬ 
sume that the residual accretion begins when deuterium 
burning ends, that is when the photospheric deuterium 
abundance drops to 1/10 of its initial value, and indicate 
this time as race. We then apply the accretion to the mod¬ 
els, starting at t=racc and assume that it declines with time 
following a power law 

M = Mof —) " [Me/yr] (3) 

where Mq is the initial rate and rj is the parameter that spec¬ 
ify the rate of decline. In principle, the parameters Mo, Tacc, 
and rj could be treated as adjustable parameters to be con¬ 
strained with the observed rates, after properly considering 
a detailed PMS evolution that includes the main accretion 
phase. However, since in this paper we perform an explo¬ 
rative analysis without a full description of the protostar 
phase and the initial evolution with large accretion rates, 
we assume the residual accretion begins when deuterium 
burning ends as defined before. 

Table [T] shows the age at the beginning of the resid¬ 
ual accretion (race) for different stellar masses. We note 
that the duration of the early large accretion, before race, 
is around 1 — 2 x 10® yr, hence negligible compared to 
the PMS lifetime (several 10^ yr). As to the initial accre¬ 
tion rate Mo, we assume Mo = 2 x 10“® Mq/jt, which 
is a reasonable value clo se to the observed m edian accre¬ 
tion rate for young TTS dHartmann et al]ll998l l. The expo¬ 


nent rj describes how fast the accretion rate declines with 
increasing age t. Her e we set r] — 0.9 to rec o ver most of the 
observ e d rates from lEsp aillat e t al.l (120141') : iGallardo et al.l 
(l2012tl : lDe March! et al) (l201lh . 


We also assume that during this residual accretion 
phase the original material is accreted following Equation 
[21 irrespective of the geometry of the disk. 

The right panel of Fig. [2] illustrates the effect of late 
accretion on the Li evolution, for a star with initial mass 
mo = O.75M0. The accretion material falling onto the star 
contains lithium with the initial abundance and, even if ac¬ 
cretion is very small after tens of million years, it restores 
the surface ^Li towards the initial value. 


2.3 EUV photo-evaporation 


Late accretion will last until the remaining gas reservoir is 
consumed or until some feedback mechanism from the star 
itself is able to clean the nearby disk. In this respect. Ex¬ 
tremely UV (EUV) radiation (13.6 — 100 eV, 10 — 121 nm) 
photons could be particular important in determining the 
end of the accretion process because they have energy high 
enough to heat the disk surface gas. The warm gas can es¬ 
cape from the gravita tional potential of the c entral star and 
flow away in a wind ijPullemond et aDl2007tl . Although no 
actual evaporative flow observation is available yet, there 
Is a large consensus that this radiation could significantly 
reduce the disk mass and shorten the lifetime of the disk. 
The existence of pre-transitional and transitional disks with 
a gap or a ho le in the disk are taken as possible evidence o f 
this process (iDullemond et aP l2007l : lEsDaillat et al.l l2014ll . 
The magnitude of the disk mass loss caus ed by EUV evap¬ 
orati on is given by the following relation llDullemond et aP 
1200711 : 


Meuv 




[Moyr (4) 


where ^euv is the EUV photon luminosity [photons/s] pro¬ 
duced by the central star. This effect can be easily included 
in PARSEC under the assumption that the stars emit as a 
black body at the given effective temperature Teff In Fig. [4] 
we show the EUV luminosity evolution for stellar models 
with initial masses mo = 0.55,0.65,0.75, and 0.85 Mq. The 
photon luminosity ^euv varies with the effective tempera¬ 
ture and the radius of the central star, therefore being im¬ 
portantly affected by the stellar evolutionary phase. Most 
of the EUV luminosities range from ^euv ~ 10®® pho¬ 
tons/s to ~ 10^^ photons/s. The residual accretion does not 
change these values much because the residual accretion rate 
is never high enough to significantly increase the mass (and 
hence the evolution) of the star. 

The EUV photo-evaporation rates obtained from 
Eq. (m are shown in Fig. [S] (dot-dashed lines) for a grid 
of low-mass stars, and they are compared with the corre¬ 
sponding accretion rates (solid lines) obtained from Eq. (|3)). 
The EUV photo-evaporation rates increase with increasing 
initial mass. Furthermore, at any given initial mass they de¬ 
pend mainly on the evolution of the EUV luminosity because 
the total mass of the star is not significantly affected by the 
accretion rate. At relatively larger masses there is an mini¬ 
mum caused by the initial shrinking of the radius at almost 
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Figure 4. Evolution of the EUV photon luminosity from the 
central stars (color-coded according to logarithmic scale shown 
on the top bar) as a function of the effective temperature and 
radius, assuming black-body emission. Note the different units 
for the radius on the left and right Y-axes. 
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Figure 5. Evolution of the residual accretion rate during the 
first 100 Myr along the PMS evolution (solid lines) Different colors 
refer to different initial masses (mo) as shown in the top color-bar. 
The residual accretion begins with the rate 2 x 10“® M^/yv and 
drops following the Eq. 0. The EUV photo-evaporation (dash- 
dotted line) rates, calculated with Eq. 0 are also plotted. From 
bottom to top are stars with initial mass 0.50Mq to 0.85 Mq as 
labeled in the color-bar. The asterisks mark the times when the 
accretion is terminated by the evaporation. 



constant effective temperature (Fig. |4]). Thereafter, the star 
evolves at constant luminosity toward the ZAMS and, since 
the effective temperature increases, the fractional number of 
EUV photons increases rapidly, so does the evaporation rate. 
When the evaporation rate is larger than the accretion rate 
the accretion effectively stops. This point is marked with an 
asterisk in the figure and the terminal age is shown in ta¬ 
ble [T] Beyond this point the star evolves at constant mass. 
Figure [5] is meant to provide a schematic diagram of this 
mechanism for the first 100 Myr. 


We would like to emphasize that the quenching of accre¬ 
tion is expected to be much more complex than our model 
description. For instance, we notice that only a fraction of 
the total EUV photons from the central star could reach 
the residual disk, if the latter maintains a small geometrical 
cross section. On the other hand, stellar activity in young 
star (e.g., accretion shock, magnetic field driven chromo¬ 
sphere activities, etc.), which is not included In our model, 
could also be a source of EUV photons and contribute to 
evaporation. This effect could balance the geometrical loss of 
EUV photons though the relative contribution between the 
two sources is not clear. The end of the late accretion phase 
could also depend on the mass of the residual disk before the 
EUV evaporation mechanism becomes effective. For metal- 
rich star-forming regions, disk masses between 0.01 — 0.2 Mq 
during the burs t accre tion pha se have been e s timat ed by 
iHartmann et al.l ll 19981 '). While lBodenheimer.1 (1201 il l esti¬ 
mate a disk masses from 0.5 Mq down to 0.0001 Mq in 
solar mass YSOs of the Taurus and Ophiuchus star-forming 
regions. The residual disk masses we are considering for the 
present work are even more uncertain than those of the main 
disk. Here we adopt the working hypothesis that the disk 
continues to provide material to the central star until EUV 
evaporation stops it. 


2.4 A combined PMS model 

A model that combines all the three effects previously dis¬ 
cussed is shown in Fig. [3] The left panel shows the PMS 
evolutionary track of a star with mass mo = 0.75Mq in the 
HR diagram. The dotted line showmass mo = 0.75M© in 
the HR diagram. The dotted line shows the PMS evolution 
at constant mass ups the PMS evolution at constant mass 
up to deuterium exhaustion. The latter point represents the 
end of the stellar birth-line, when the phase of large accre¬ 
tion terminates and the star begins the contraction at almost 
constant mass downward in the HR diagram. The solid line 
begins when the deuterium abundance reduces to 1/10 of 
its initial value. At this time we apply an accretion with a 
rate starting at 2 x 10“® MQ/yr (dark blue line in the right 
panel). This value is too low to prevent a rapid contraction 
and the central temperature Teen (brown dashed line) rises 
until surface Li depletion begins (red line). 

At this stage the star is fully convective and the sur¬ 
face Li depletion proceeds faster than the restoring effect of 
the accretion because of the high central temperature. After 
about five Myr, the core becomes radiative and the tem¬ 
perature at the base of the convective region Ti,cz (brown 
solid line) begins to decrease followed by a decline of the Li 
burning rate. Then the photospheric Li abundance begins 
to increase again, especially when Th^z falls below 4 x 10® K. 
This is because on one side Li nuclear burning quenches off 
and, on the other, the dilution of the infalling Li becomes 
weaker as the mass contained in the the convective envelope 
becomes smaller. The mass accreted in this phase is rela¬ 
tively little and only slightly affects the total mass of the 
star with negligible effects on Its structure. In practice, the 
star evolves at constant mass. 

The surface Li abundance steadily increases towards its 
primordial value, and will reach It If accretion does not ter¬ 
minate. This is due to the shrinking of the surface convective 
zone which becomes so small that even a residual accretion 


© 0000 RAS, MNRAS 000, 000-000 




























6 X. Fu, et al 


Table 1. Relevant parameters for stars with different initial masses. Tacc is the age at the beginning of the late accretion phase; m, is 
stellar mass in main sequence; tend is the age at the end of residual accretion; An{Li) p mS and An{Li) mS denote the decrement in the 
Li number density during pre-main sequence and main sequence, respectively. Note that the PMS Li depletion is stronger than the MS 
one. 


mo {Mq) 

0.57 

0.58 

0.59 

0.60 

0.62 

0.65 

0.70 

0.75 

0.80 

Tacc (10®J/r) 

1.876 

1.873 

1.834 

1.787 

1.731 

1.689 

1.595 

1.560 

1.433 

m* {Mq) 

0.616 

0.625 

0.633 

0.641 

0.658 

0.685 

0.729 

0.775 

0.821 

tend {Myr) 

585.55 

518.96 

460.83 

399.76 

309.84 

211.46 

110.98 

61.14 

35.46 

An{Li)pMS (10“^®) 

3.16 

2.81 

2.60 

2.50 

2.33 

2.11 

1.75 

1.17 

1.21 

An{Li)MS (10“i®) 

0.70 

0.78 

0.65 

0.49 

0.41 

0.50 

0.86 

1.19 

1.11 



Figure 6. H-R diagram from PMS with overshoot, accretion, and 
photo-evaporation, to the end of the main sequence. The tracks 
are for stars with initial mass 0.85, 0.80, 0.75, 0.70, 0.65, 0.62, 
0.60, 0.59, 0.58, 0.57, 0.55, and 0.50 Mq from top to bottom. 
Black dashed lines at the upper-right side are Hayashi line with 
constant initial masses, which is not the real case for the late 
accretion. Solid lines are the stellar evolutionary tracks after deu¬ 
terium burning with an initial lithium value A(Li)=2.72. Tracks 
are color-coded according to their surface lithium abundance, 
as indicated in the right vertical bar, starting from A(Li)=2.72 
(in red) to A (Li) = 0 (violet). The black color corresponds to 
A(Li)< 0. The area between the two vertical dotted lines is the 
effective temperature range of the Spite plateau. 


rate less than Mq/jx is enough to engulf the thin sur¬ 

face convective layers. The interplay between the size of the 
convective regions and the evolution of the accretion rates 
are critical for reproducing a given pattern in the surface Li 
abundance as a function of the initial mass on the main se¬ 
quence. At an age of about 60 Myr (see TableO, the residual 
accretion is terminated by the EUV photo-evaporation. 

The H-R diagram of the full sample of stars consid¬ 
ered in this paper is shown Fig. [6] The area between the 
two vertical dotted lines illustrates the temperature range 
of the Spite plateau. The high end of this range is also the 
temperature threshold of the observable lithium. For stars 
with higher Tefj (the more massive ones), lithium is almost 
fully ionized in the stellar photosphere. Likewise in Fig. |31 
the upper Hayashi lines before the completion of deuterium 
burning are drawn with dotted lines whilst the solid lines 


correspond to the stellar evolutionary tracks after deuterium 
burning. The surface abundance of Li is color-coded, start¬ 
ing from an initial abundance A(Li)=2.72. After the end of 
the accretion (marked with the asterisks), stars evolve with 
constant masses along main sequence and Li abundance no 
longer increases. 


3 MAIN SEQUENCE LI EVOLUTION 

During the main sequence phase Li is depleted by burning at 
the base of the convective zone because, even a low nuclear 
reaction rate at Tt,cz ~ 2 x 10® K, could cause significant 
depletion in a timescale of several Gys. However, for masses 
larger than mo = O.6OM0, Li burning is insignificant. 

Another effect that could modify the photospheric Li 
abundance in low mass POP H stars is microscopic diffu¬ 
sion. This long-term stellar process cannot be observed di¬ 
rectly in the star, and its efficiency needs to be calibrated. 
Microscopic diffusion shortens the main sequence lifetime 
and leads to a depletion of the surface elements. The solar 
model requires the inclusion of microscopic diffusion, oth¬ 
erwise even the age of the Sun could not be correctly pro¬ 
duced. We calculated evolutionary models up to the end of 
the main sequence including pressure diffusion, temperat ure 
diffusion, and concentration diffusion dThoul et alJll994l L It 
has already been shown that microscopic diffusion without 
any correction for radiative levitation and surface convec¬ 
tive turbule nce cannot reproduce the observed photospheric 
abundances llRichard et al.ll2002l ). In the case of Li, the sur¬ 
face convective zone of the stars at the plateau high temper¬ 
ature end (the relatively more massive stars) is so thin that, 
if microscopic diffusion is at work and not balanced by the 
radiative levitation, a too strong de pletion is produc e d. For 
example, according to the model of ISalaris fc Wei^ (1200 ll ) 
we can see that even at Tefi ~ 6200 K the predicted Li de¬ 
pletion is too strong compared to the observed plateau. 

This does not happen with the standard PARSEC code 
because gravitational settling is always inhibited in the out- 
ermos t region of the enve lope, AM — 0.5% of the stellar 
mass llBressan et al.ir2013l ). T his choice of AM is m ade ac¬ 
cording to the suggestion by IChabover et ^ ll200ll) . They 
noticed that, while observations indicate that the relatively 
hot low mass stars at the turnoff of the globular cluster 
NGC 6397 show th e same surface [Fe/H ] abundance of the 
evolved RGB stars (iGratton et al.l[200ll ). models with grav¬ 
itational settling that well reproduce the solar data predict 
that they should show a surface [Fe/H] abundance at least 
0.28 dex lower. They thus conclude that gravitational set- 
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Figure 7. A(Li) as a function of the stellar age for stars with different initial masses. Different colors refer to different initial masses 
(mo). The horizontal grey dashed line is the SBBN prediction (A(Li)=2.72), while the horizontal grey solid one indicates A(Li) for the 
Spite plateau. Left panel: Li abundance during the PMS phase during the first 100 Myr. The initial mass (mo) of the models decreases 
from 0.85 Mq to 0.50 Mq, from top to bottom in the rapid rising branch on the left of the diagram, as indicated in the legend. The 
asterisks mark the end of the accretion. Right panel: Li evolution up to an age of 12 Gyr, zooming in around the region around the Spite 
plateau. Color codes and the line styles are the same as in the left panel. Note that stars with initial masses greater than 0.80 Mq, have 
a main sequence life shorter than 10 Gyr, and they may have already evolved. 


tling is inhibited in the outermost layers of such stars and 
show that, in order to reconcile the predicted with the ob¬ 
served [Fe/H] abundance, the size of this layer should be 
A M ~ 0.5% — 1% Mq for a star of 1 Mq. We checked 
that changing the standard PARSEC parameter in the range 
from AM — 0.1% to A M = 1% does not appreciably affect 
our results. 


4 RESULTS 

We analyse the importance of overshoot and late accre¬ 
tion during the pre-main sequence phase of stellar evolution 
by computing a set of evolutionary tracks of initial stellar 
masses 0.85, 0.80, 0.75, 0.70, 0.65, 0.62, 0.60, 0.59, 0.58, 0.57, 
0.55, and 0.50 Mq, assuming an efficient overshoot, as de¬ 
scribed in Sect. and exploring different values for the 
parameters that describe the form of the residual accretion 
rate, Mq and rj. We discuss below the results obtained as¬ 
suming Mo = 2 X lO“®M0/yr and rj = 0.9 in equation |31 
which produces accretion rates compatible with the observed 
values. The models are calculated from the PMS till the end 
of the main sequence. 

We will first focus on the evolution of Li abundance 
during the PMS phase. The surface Li abundance evolution 
of our selected models during the first 100 Myr of PMS phase 
is plotted in the left panel of Fig. [T] With our choice of the 
parameters, the initial Li evolution in these stars is regulated 
by overshooting. The convective overshooting is so efficient 
that the photospheric Li is significantly depleted after the 
first few Myr. The depletion caused by overshooting gets 


stronger at the lower masses. At later times the core becomes 
radiative and the convective envelope starts receding. 

At subsequent stages Li evolution is governed by the 
competition between late accretion, dilution and nuclear 
burning at the bottom of the convective envelope. As soon 
as efficient Li burning ceases, accretion begins to restore the 
surface Li abundance. In the more massive stars the con¬ 
vective zone shrinks more rapidly and Li restoring is faster. 
Conversely, in lower mass stars where Li depletion is more ef¬ 
ficient because of the deeper convective zones, the Li restor¬ 
ing time is longer. 

Without EUV evaporation, more massive stars rapidly 
recover the initial surface Li abundance while in lower mass 
stars this process is slower because nuclear reactions will 
continue to burn the accreted Li and the convective dilution 
is larger. The effects of EUV evaporation strongly modify 
this picture and is regulated by the stellar mass. More mas¬ 
sive stars evolve more rapidly toward the ZAMS and their 
EUV luminosity is able to stop accretion at earlier times, 
even before the initial abundance is completely restored (see 
Fig. 0. On the other hand, in lower mass stars where Li 
depletion is more pronounced, the EUV evaporation rates 
are lower and accretion, though decreasing, is still able to 
drive the surface Li abundance toward the initial value, un¬ 
til It is inhibited by evaporation. At even lower masses (e.g., 
mo = 0.50 Mq), the recovery of the initial Li abundance is 
delayed or even inhibited by the deep convective zone which 
is still large when the model reaches the ZAMS. In this case 
Li restoring is prevented both by an efficient dilution, and 
possibly by the nuclear burning at the base of such deep 
convective envelopes. 
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Figure 8. Our results in comparison with the lithium abundance 
measuremen ts in POP II stars. T he grey filled squares are POP 
II data from l^olaro et al.l ll2012h . Our predictions are shown for 
stars at the end of the late accretion phase (filled triangles), and 
on the main sequence at 10 Gyr (filled upper circle) and 12 Gyr 
(filled lower circle). Symbols are color-coded according to the ini¬ 
tial stellar mass, from the left to the right are 0.80 Mq to 0.57 
Mq as the same as labeled. The black dashed line marks the 
primordial Li abundance according to the SBBN. 

In summary, our analysis indicates that Li evolution 
during the pre-main sequence phase may be regulated by 
the competition between the destroying effect of overshoot¬ 
ing at early times, and the restoring effect of residual accre¬ 
tion, which is in turn modulated by EUV photo-evaporation 
quenching. This latter effect controls the end result, depend¬ 
ing on the mass. At the higher stellar masses Li restor¬ 
ing is faster, but accretion quenching takes place earlier. 
At lower stellar masses the time-scale for Li restoring is 
longer, but the accretion lasts longer because the EUV 
photo-evaporation has a lower efhciency. The combination 
of these effects in stars with different initial masses tends to 
level out the Li abundance just after the PMS at a value 
that is already below the initial one. 

Galactic h alo POP II stars we observe today are about 
10-12 Gyr old (Ijofre fc WeissiBoill ). In order to reproduce 
the observed Li abundance of these stars we have evolved 
our models up to the end of central hydrogen burning (till 
the turn-off phase). In this phase we also account for micro¬ 
scopic diffusion which is known to be a long term effect that 
can modify the photospheric element abundances. The evo¬ 
lution of the photospheric Li abundance during the main se¬ 
quence is depicted in the right panel of Fig. [T] The predicted 
photospheric Li abundances are compared with those of the 
Galactic halo POP II stars in Fig.[Sl We plot the abundances 
at the end of the accretion phase (filled triangles) and those 
at older ages, 10 Gyr (filled upper circles) and 12 Gyr (filled 
lower circles). Stars with initial mass mo ^ 0.85Mq are ex¬ 
cluded because their effective temperatures on the MS are 
warmer than the observed ones (see Fig. and also because 
their main sequence lifetimes are shorter than the relevant 
age range. For our standard choice of parameters, stars with 
initial mass from mo = 0.62Mq to O.SOM©, nicely populate 
the Spite plateau (A(Li) ~ 2.26). 

Models with lower mass, mo = 0.57Mq to 0.60M©, fall 
on the observed declining branch towards lower tempera- 


Figure 9. By applying the same parameter of envelope over¬ 
shooting and accretion, we could reproduce the Spite plateau and 
the first Li decline branch for a wide range of metallicities (from 
Z=0.00001 to Z= 0.0005). The compared POP II data are the 
same as in figure |8] The model results are all main sequence stars 
at age 10 Gyr, with initial mass 0.80 M© to 0.57 M© from the left 
to the right for each metallicity as shown in the color-bar label. 

tures. We confirm that the first part of this branch is pop¬ 
ulated by low-mass main sequence stars. Indeed this is the 
signature of the strong Li depletion during the PMS phase 
(filled triangles), followed by further depletion during the 
main sequence evolution (filled upper and lower circles). 

Both PMS Li depletion and MS diffusion contribute to 
make the total A(Li) decrease, with the former process play¬ 
ing the main role. Table [T] lists the Li number density de¬ 
crease during PMS and MS, the former is much more than, 
or at least the same as, the latter, especially for the lower 
masses. The first part of Li declining branch is caused mainly 
by the PMS depletion. 

The Spite plateau is populated by metal-poor main se¬ 
quence stars with different metallicities. With the same over¬ 
shooting and accretion parameters we can indeed reproduce 
the plateau, including the initial Li declining branch, over a 
wide range of metallicities, (from Z=0.00001 to Z=0.0005, 
that is from [M/H]=-3.2 to [M/H]=-1.5) as shown in Fig. 


5 DISCUSSION AND CONCLUSION 

In the previous sections we have shown how the effects of 
efficient envelope overshooting, residual accretion during the 
pre-main sequence phase, and microscopic diffusion during 
the main sequence may modify the photospheric Li abun¬ 
dance in low-metallicity stars. We have also considered the 
EUV photo-evaporation process which, by terminating the 
accretion phase, could introduce a sort of self-regulating pro¬ 
cess. 

In order to give a more exhaustive picture of the im¬ 
portance of the various processes we computed more sets 
of models under different assumptions concerning the above 
physical mechanisms. We synthesize the results of these ad¬ 
ditional models In Fig. 1101 where we compare the predicted 
Li abundances with the observed ones (see Fig. EJ. 

• Standard models 

We begin with two sets of models where none of the PMS 
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Figure 10. Comparison of the standard models and the present 
mo dels for metallici t y Z= 0.0001 together with the POP II data 
and lSalaris &: Weisj i200lll pure diffusion model. Different sets of 
models are classified in the legend. The black d ash dot dot line is 
extracted from Fig. 1 of ISalaAs &: Wei^ ll200lh with [Fe/H]=-2.6 
and age 12 Gyr. 


effects have been considered. The first set (orange dashed 
fine) is the standard model without microscopic diffusion 
while, In the second set of model (green dot-dashed line) 
microscopic diffusion is introduced as it is in the standard 
PARSEC model. 

The model without diffusion is clearly at variance with 
the observed data. On the other hand when diffusion is in¬ 
cluded the Li abundance reaches A(Li)=2.45 dex, produc¬ 
ing a plateau somewhat higher than the o bserved one, by 
a fact or of 1.6. The pure diffusion model of ISalaris fc Wei^ 
ll 200 lll is also plotted in this figure (extracted from figure 1 
of Salaris fc Weiss! ll200lll with [Fe/H]=-2.6 and age 12 Gyr). 
Contrary to them, our diffusive model does not produce the 
strong downturn at high effective temperatures. This is due 
to the inhibition of diffusion in the outer layers of the star 
which, as discussed in Sect. |31 is required to prevent a too 
strong sedimen tation of helium and heavy elements at the 
stellar surface llRichard et al.|[2002l : iBressan et al.ll2012l L 

While the plateau could eventually be reproduced by in¬ 
creasing the efficiency of diffusion strongly, we note however 
that, with the adopted d iffusion parameters the solar model 
is very well reproduced (IBressan et al.ll2012f ). 

• Present models 

The present models, that include PMS effects and diffu¬ 
sion (red solid fines in Fig. llOH . provide a good fit to the data. 
They reproduce both the observed plateau and the declin¬ 
ing branch at low temperatures. As a test, we recomputed 
the same models but switching off microscopic diffusion dur¬ 
ing the main sequence phase (blue dotted line). These new 
models are identical in the PMS evolution and incorporate 
the Li nuclear burning during main sequence. The latter is 
mainly responsible for the behavior of Li depletion at low 
effective temperatures. They do not reach the plateau value, 
terminating somewhat at higher Li abundances, suggesting 
that some additional depletion is required. 


In our model we do not consider stellar rotation. The 


additional effects of rotation on Li depletion have been dis¬ 
cussed in many works, but mainly in the context of solar- 
type stars. During the PMS, rotation may affect the internal 
structure of the stars both because of the presence of the 
centrifugal force and the effects of rotational mixing at the 
base of the convective envelope. lEggenberger et aU ll2012h 
conclude that by including the full treatment of rotation 
in PMS, 0.1 dex more Li depletion is obtained. During the 
evolution on the main sequence, rotational mixing tends to 
oppose to gravitational settling, and this solo effect should 
decrease the surface elements depletion with respect to non¬ 
rotating models. However in the case of Li, rotational mixing 
may be so efficient that Li is dragged to encounter temper¬ 
atures which could efficiently destroy it by nuclear burning. 
Th is is the case of the slow ro tator model of M= 1. 0 Mp , 
of lEggenberger et ahl (|201(]ll . iGarcIa Lopez et aP (Il994fi 
find that because of the rotation-caused structural differ¬ 
ence, stars with faster rotation in open clusters Pleiades 
and a Persei have higher Li abundance than the slower 
ro taters. This r esults is confirmed by Pleiades observation 
of I Jones et al.l (Il997f l. In general it is found that rotation 
leads to a large spread in the main sequence Li abundance 
of solar-type stars a t varying mass jMartin fc Clareull996l : 
iMendes et al.l 119991 1 , which could be difficult to reconcile 
with that observed in the Spite plateau of metal-poor stars. 

It is important to stress that, as long as the stars are 
accreting at significant rates, deuterium burning keep their 
central temperature around 10®K. During this phase the 
stars evolve along the birth line increasing their masses (e.g. 
reaching 1 Mq in 0.1 Myr for an accretion rate of 10~® 
Mq /yr as indi cated in table [H), while preserving the initial 
Li abundance flStahleilll983ll . Once the accretion rates fall 
below several 10 “® M© /yr the stars abandon the stellar 
birth-line and evolve at almost constant mass. It is from 
this point that we consider the additional effects of envelope 
overshooting, which reduces the Li abundance, and of accre¬ 
tion, which tends to restore the original abundance. We have 
also shown that EUV evaporation of material falling into the 
stars can act as a natural self-regulating mechanism. 

In the relatively more massive stars, where the original 
Li abundance could be more quickly restored, EUV evapora¬ 
tion is more efficient and quenches out accretion earlier than 
in lower mass stars, which thus have more time to accrete the 
pristine gas. These combined effects tend to level out the Li 
abundance in stars of different masses. At even lower masses, 
Li restoring is inhibited by nuclear burning in the deep con¬ 
vective envelopes. After the PMS, Li is slowly depleted by 
microscopic diffusion and nuclear burning which act during 
main sequence phase. For stars with mo 15 O.62M0 diffu¬ 
sion drives the main Li depletion whilst, for mo ^ O. 6 OM 0 
stars, LI is significant burned at the base of the convective 
envelope. 

In this scenario, the observed Spite plateau and its 
falling branch at low temperatures are well reproduced by 
our models, within a reasonable parametrization of the late 
accretion mechanism. Very efficient overshoot is only ap¬ 
plied during the early PMS to favour strong Li depletion at 
the beginning of the evolution. We also find that, even for 
a maximum envelope overshoot Ae = 0.7 Hj,, which is th e 
value calibrated from the RGB bump (lAlongi et al.|[l99lll . 
the plateau and its falling branch can be finally well re¬ 
produced for the same accretion rate. This indicates that 
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a wide range of envelope overshoot could be applied to the 
PMS phase. Even for Ae = 0.7 Hp the PMS Li depletion 
is very strong; without accretion a very low Li abundance 
would be observed, e.g. for mo = O.IOMq A(Li) drops to 1 
dex before Li restoring. It would be important to directly 
observe whether the very early Li depletion phase exists 
or not, because it could confirm that efficient overshoot is 
at work. Indeed, at higher metallicities there is evidence 
of a si gnificant drop of the Li abun dance during the PMS 
phase dSomers &: Pinsonneaultll2014ll . On the contrary, the 
lack of nearby metal-poor star forming regions prevents any 
observational test at low metallicity. This opportunity will 
be perhaps offered by the next generation telescopes (E- 
ELT, TMT, and GMT) which could observe the details of 
low mass metal-poor PMS stars in star forming regions like 
those of the S agittarius dwarf irreg ular galaxy (metallicity 
at Z ~ 0.0004 iMomanv et al.l[2002 fL 

For stars with initial mass mo = 0.62 Mq — 0.80 Mq, 
ages in the range 10-12 Gyr, and a wide range of metallicities 
(Z=0.00001 to Z=0.0005), our final model is able to predict 
a present day abundance consistent with the Spite plateau, 
although starting from an initial value of A(Li)=2.72 as 
inferred from the baryon-to-photon ratio suggested by the 
GMB and by the deuterium measurements. Thus it offers 
a likely mechanism to solve the long standing lithium prob¬ 
lem. It also reproduces the observed lithium drop at the low- 
temperatures. The proposed solution relies on stellar physics 
and evolution, assuming the validity of the current SBBN 
theory. Though this model is quite schematic and rough, 
it clearly suggests that the Spite plateau could be the net 
result of mechanisms more complex than those considered 
so far. In particular, a key role may be attributed to the 
interaction between the stars and their environments. 

If the lithium abundance we presently observe in POP II 
stars has been restored in the stellar atmospheres by a resid¬ 
ual tail of accretion after a phase of strong depletion, then 
other puzzling observations could be explained, namely: 

• Li at very low meta llicities. 

ISbordone et al.l (l2010l ') found in their observations that 
for the extreme m etal-poor sta r s, Li abun dance drops 
when [Fe/H]< —3. iGaffau et al.l ll201ll ') and iFrebel et al.l 
(l2005ll observed dwarf s with [Fe/H] at the l owest levels with¬ 
out any detectable Li. iHansen et al. lim detected Li in an 
un-evolved star with [Fe/H]=-4.8 at a level of A(Li) = 1.77 
and anther recently discovered extremely metal poor dwarf 
SDSS J1742+2531 with [Fe/H] = -4.8 shows A(Li) < 1.8 
llBonifacio et al.ll2015l L The low Li abundance in these stars 
could be explained by a failed or weaker late accretion. For 
the most metal-poor stars, late accretion during the PMS 
phase might be inhibited and Li could not be restored. As 
we have already shown in Fig. [2] (right panel), if no accretion 
takes place after the initial Li is depleted by the convective 
overshoot, Li in these stars would correspond to a very low, 
or even undetectable, abundance. Thus the model presented 
here provides a possible key to account for this phenomenon, 
and we will discuss it in detail in a following paper. 

• Spectroscopic binaries and outliers 

Our model might also provide an explanation to 
the Li discrepancy observed in metal-poor binary 
stars. The POP II spectroscop ic binaries GS 22876- 
032 li Gonzales Hernandez et al.f l2008h and G 166- 


45 llAoki et al.l I 2 OI 2 II are composed by dwarfs with 
effective temperatures characteristic of the Spite plateau, 
but with the primary stars showing slightly higher Li 
abundances. It could simply be an effect of the competition 
between the components in the late PMS accretion, with 
the primary star being the favorite one. 

The failure or increase of the accretion process could also 
provide an explanation for the few POP H ^Li depleted stars 
or the few POP H stars that exhibit a ^Li abundances at the 
SBBN prediction level, respectively. Among the latter there 
is BD +2 3 3912 with A(Li)=2 .60 which stands out from the 
others in I Bonifacio fc Molard lll997tl . 

• Lithium behavior in POP I stars 

Li abundance at the solar system for mation time, as ob¬ 
taine d from meteorites, is A(Li)= 3.34 llAnders &: Grevessd 
Il989ll . The hot F sta rs of young open clusters never reach 
this meteoritic value (iFord et al. 1200 il l despite they should 
have an initial Li value even higher than the meteoritic one 
due to the Galactic Li increase with time. This Li drop can¬ 
not be explained by main sequence Li evolution because the 
stars are very young. A PMS Li modihcation could be the 
possible mechanism for the depletion. This will be discussed 
in detail in a following work dedicated to the metal-rich 
stars. 

• Low lithium abund anc e in the planet h ost stars 

llsraelian et al.l (l2009h and iGonzaled ll2014h present the ob¬ 
servational evidence that the solar analogues with planet (s) 
show lower lithium abundanc e than those without dete cted 
planets. To explain this effect f^genberger et al.l 1 I 2 OIOII ad¬ 
vance the hypothesis that giant exoplanets prefer stars with 
slow rotation rates on the ZAMS due to their longer disk life¬ 
time during the PMS evolution. Then their ZAMS models 
with rotational mixing, which we already discussed in sec¬ 
tion |31 would indicate an enhanced lithium depletion, con¬ 
sistent with observations. Again it is worth noting that a 
strong PMS depletion at decreasing mass is predicted b y 
solar metallicity models with rotation (iMendes et al.|[T999ll . 
and it is not ye t clear how this could be r econciled with th e 
observations of llsraelian et all ll2009ll and iGonzale j ll2014ll . 
In our model, this behavior could be explained if somehow 
the planet (s), located within the accreting disk, interfere 
with the late accretion process. The PMS accretion could 
end before being terminated by the EUV evaporation, in¬ 
hibiting efficient lithium restoring. Thus these stars will 
show lower Li abundance than their analogues without de¬ 
tected planets. This effect could also explain the larger Li 
abundance scatter observed in metal-rich st ars, which are 
more likely to host planets (llda fc Linll2004l L with respect 
to what observed in metal poor stars. 
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